Additive domain decomposition operator splittings -- convergence analyses in a dissipative framework

We analyze temporal approximation schemes based on overlapping domain decompositions. As such schemes enable computations on parallel and distributed hardware, they are commonly used when integrating large-scale parabolic systems. Our analysis is conducted by first casting the domain decomposition procedure into a variational framework based on weighted Sobolev spaces. The time integration of a parabolic system can then be interpreted as an operator splitting scheme applied to an abstract evolution equation governed by a maximal dissipative vector field. By utilizing this abstract setting, we derive an optimal temporal error analysis for the two most common choices of domain decomposition based integrators. Namely, alternating direction implicit schemes and additive splitting schemes of first and second order. For the standard first-order additive splitting scheme we also extend the error analysis to semilinear evolution equations, which may only have mild solutions.Comment: Please refer to the published article for the final version which also contains numerical experiments. Version 3 and 4: Only comments added. Version 2, page 2: Clarified statement on stability issues for ADI schemes with more than two operator

A full space-time convergence order analysis of operator splittings for linear dissipative evolution equations

The Douglas--Rachford and Peaceman--Rachford splitting methods are common choices for temporal discretizations of evolution equations. In this paper we combine these methods with spatial discretizations fulfilling some easily verifiable criteria. In the setting of linear dissipative evolution equations we prove optimal convergence orders, simultaneously in time and space. We apply our abstract results to dimension splitting of a 2D diffusion problem, where a finite element method is used for spatial discretization. To conclude, the convergence results are illustrated with numerical experiments

Spatial and Physical Splittings of Semilinear Parabolic Problems

Splitting methods are widely used temporal approximation schemes for parabolic partial differential equations (PDEs). These schemes may be very efficient when a problem can be naturally decomposed into multiple parts. In this thesis, splitting methods are analysed when applied to spatial splittings (partitions of the computational domain) and physical splittings (separations of physical processes) of semilinear parabolic problems. The thesis is organized into three major themes: optimal convergence order analysis, spatial splittings and a physical splitting application.In view of the first theme, temporal semi-discretizations based on splitting methods are considered. An analysis is performed which yields convergence without order under weak regularity assumptions on the solution, and convergence orders ranging up to classical for progressively more regular solutions. The analysis is performed in the framework of maximal dissipative operators, which includes a large number of parabolic problems. The temporal results are also combined with convergence studies of spatial discretizations to prove simultaneous space–time convergence orders for full discretizations.For the second theme, two spatial splitting formulations are considered. For dimension splittings each part of the formulation represents the evolution in one spatial dimension only. Thereby, multidimensional problems can be reduced to families of one-dimensional problems. For domain decomposition splittings each part represents a problem on only a smaller subdomain of the full domain of the PDE. The results of the first theme are applied to prove optimal convergence orders for splitting schemes used in conjunction with these two splitting formulations. The last theme concerns the evaluation of a physical splitting procedure in an interdisciplinary application. A model for axonal growth out of nerve cells is considered. This model features several challenges to be addressed by a successful numerical method. It consists of a linear PDE coupled to nonlinear ordinary differential equations via a moving boundary, which is part of the solution. The biological model parameters imply a wide range of scales, both in time and space. Based on a physical splitting, a tailored scheme for this model is constructed. Its robustness and efficiency are then verified by numerical experiments

Efficient simulations of tubulin-driven axonal growth

This work concerns efficient and reliable numerical simulations of the dynamic behaviour of a moving-boundary model for tubulin-driven axonal growth. The model is nonlinear and consists of a coupled set of a partial differential equation (PDE) and two ordinary differential equations. The PDE is defined on a computational domain with a moving boundary, which is part of the solution. Numerical simulations based on standard explicit time-stepping methods are too time consuming due to the small time steps required for numerical stability. On the other hand standard implicit schemes are too complex due to the nonlinear equations that needs to be solved in each step. Instead, we propose to use the Peaceman--Rachford splitting scheme combined with temporal and spatial scalings of the model. Simulations based on this scheme have shown to be efficient, accurate, and reliable which makes it possible to evaluate the model, e.g.\ its dependency on biological and physical model parameters. These evaluations show among other things that the initial axon growth is very fast, that the active transport is the dominant reason over diffusion for the growth velocity, and that the polymerization rate in the growth cone does not affect the final axon length.Comment: Authors' accepted version, (post refereeing). The final publication (in Journal of Computational Neuroscience) is available at Springer via http://dx.doi.org/10.1007/s10827-016-0604-

Analyses and Applications of the Peaceman--Rachford and Douglas--Rachford Splitting Schemes

Splitting methods are widely used as temporal discretizations of evolution equations. Such methods usually constitute competitive choices whenever a vector field can be split into a sum of two or more parts that each generates a flow easier to compute or approximate than the flow of the sum. In the research presented in this Licentiate thesis we consider dissipative evolution equations with vector fields given by unbounded operators. Dynamical systems that fit into this framework can for example be found among Hamiltonian systems and parabolic and hyperbolic partial differential equations (PDEs). The goal of the presented research is to perform convergence analyses for the lternating direction implicit (ADI) methods in the setting of dissipative operators. In this context these methods are known to possess excellent stability properties. Additionally, they generate easily computable numerical flows and are ideal choices for studying convergence to stationary solutions, a property related to their favorable local error structure. In this thesis we consider the Peaceman--Rachford and Douglas--Rachford schemes, which were the first ADI methods to be constructed and to this day are the most representative members of the ADI method class. We perform convergence studies for the Peaceman--Rachford and Douglas--Rachford schemes when applied to semilinear, dissipative evolution equations, that is, when the vector fields are given by the sum of a linear and a nonlinear dissipative operator. Optimal convergence orders are proven when the solution is sufficiently regular. With less regularity present we are still able to prove convergence, however of suboptimal order or without order. In contrast to previous convergence order analyses we do not assume Lipschitz continuity of the nonlinear operator. In the context of linear, dissipative evolution equations we consider full space-time discretizations. We assume that the full discretization is given by combining one of the two aforementioned ADI methods with a general, converging spatial discretization method. In this setting we prove optimal, simultaneous space-time convergence orders. Advection-diffusion-reaction models, encountered in physics, chemistry, and biology are important examples of dissipative evolution equations. In this thesis we present such a model describing the growth of axons in nerve cells. The model consists of a parabolic PDE, which has a non-trivial coupling to nonlinear ordinary differential equations via a moving boundary, which is part of the solution. Since additionally the biological model parameters imply a wide range of scales, both in time and space, the application of a numerical method is involved. We make an argument for a discretization consisting of a splitting which is integrated by the Peaceman--Rachford scheme. The choice is motivate by the results of some numerical experiments

The Center for Integrated Molecular Brain Imaging (Cimbi) database

AbstractWe here describe a multimodality neuroimaging containing data from healthy volunteers and patients, acquired within the Lundbeck Foundation Center for Integrated Molecular Brain Imaging (Cimbi) in Copenhagen, Denmark. The data is of particular relevance for neurobiological research questions related to the serotonergic transmitter system with its normative data on the serotonergic subtype receptors 5-HT1A, 5-HT1B, 5-HT2A, and 5-HT4 and the 5-HT transporter (5-HTT), but can easily serve other purposes.The Cimbi database and Cimbi biobank were formally established in 2008 with the purpose to store the wealth of Cimbi-acquired data in a highly structured and standardized manner in accordance with the regulations issued by the Danish Data Protection Agency as well as to provide a quality-controlled resource for future hypothesis-generating and hypothesis-driven studies.The Cimbi database currently comprises a total of 1100 PET and 1000 structural and functional MRI scans and it holds a multitude of additional data, such as genetic and biochemical data, and scores from 17 self-reported questionnaires and from 11 neuropsychological paper/computer tests. The database associated Cimbi biobank currently contains blood and in some instances saliva samples from about 500 healthy volunteers and 300 patients with e.g., major depression, dementia, substance abuse, obesity, and impulsive aggression. Data continue to be added to the Cimbi database and biobank

Palmitate-Induced β-Cell Dysfunction Is Associated with Excessive NO Production and Is Reversed by Thiazolidinedione-Mediated Inhibition of GPR40 Transduction Mechanisms

BACKGROUND: Type 2 diabetes often displays hyperlipidemia. We examined palmitate effects on pancreatic islet function in relation to FFA receptor GPR40, NO generation, insulin release, and the PPARgamma agonistic thiazolidinedione, rosiglitazone. PRINCIPAL FINDINGS: Rosiglitazone suppressed acute palmitate-stimulated GPR40-transduced PI hydrolysis in HEK293 cells and insulin release from MIN6c cells and mouse islets. Culturing islets 24 h with palmitate at 5 mmol/l glucose induced beta-cell iNOS expression as revealed by confocal microscopy and increased the activities of ncNOS and iNOS associated with suppression of glucose-stimulated insulin response. Rosiglitazone reversed these effects. The expression of iNOS after high-glucose culturing was unaffected by rosiglitazone. Downregulation of GPR40 by antisense treatment abrogated GPR40 expression and suppressed palmitate-induced iNOS activity and insulin release. CONCLUSION: We conclude that, in addition to mediating acute FFA-stimulated insulin release, GPR40 is an important regulator of iNOS expression and dysfunctional insulin release during long-term exposure to FFA. The adverse effects of palmitate were counteracted by rosiglitazone at GPR40, suggesting that thiazolidinediones are beneficial for beta-cell function in hyperlipidemic type 2 diabetes

Högstadieelevernas matvanor i Lilla Edet : vad har hänt på elva år? Vilka variationer finns? Vad blir konsekvenserna?

Syftet med uppsatsen är att undersöka matvanorna hos Lilla Edets åttonde klassare och den påverkan de har på deras vardag. Finns det ett samband mellan frukostvanor och BMI? Tidigare forskning har visat att sunda frukostvanor motverkar övervikt. Min avsikt är att se om jag kan finna liknande samband bland Lilla Edets åttondeklassare. Jag avser även att söka efter samband mellan sunda matvanor och prestationsförmåga i skolan samt vilken kunskap respondenterna har om kost. Slutligen avser jag att studera vilka skillnader som finns mellan de båda könen. För att få en uppfattning om i vilken riktning utvecklingen är på väg har jag valt GÖTE-ÄLV studien som referensobjekt. Det är en enkätundersökning som utfördes 1996 bland regionens alla åttondeklassare. Av den anledningen har jag valt att utgå ifrån samma enkät med vissa omarbetningar. De mest intressanta resultaten av undersökningen var att medelvärdet på BMI hos respondenterna låg kvar på samma nivå som i GÖTE-ÄLV studien för elva år sedan. Att det inte fanns något samband mellan frukostvanor och BMI. Att skillnaderna mellan könen fortfarande är stora. Och slutligen att andelen som äter frukost varje dag har minskat